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Statistics of Dissipation and Enstrophy Induced by a Set of Burgers Vortices 



OO 



0^ 



> 

o 
m 
o 

00 
I 

O 

^ . 



Guowei He^^^, Shiyi Chen^'^, Robert. H. Kraichnan'^ , Raoyang Zhang^, and Ye Zhou^ 

^Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, NM 87545 
IBM Research Division, T. J. Watson Research Center, P.O. Box 218, Yorktown Heights, NY 10598 
^396 Montezuma 108, Santa Fe, New Mexico, 87501-2626, 
* LNM, Institute of Mechanics, Chinese Academy of Sciences, Beijing, 100080, P. R. China 

Dissipation and enstropy statistics are calculated for an ensemble of modified Burgers vortices in 
equilibrium under uniform straining. Difi'erent best-fit, finite-range scaling exponents are found for 
locally-averaged dissipation and enstrophy, in agreement with existing numerical simulations and 
experiments. However, the ratios of dissipation and enstropy moments supported by axisymmetric 
vortices of any profile are finite. Therefore the asymptotic scaling exponents for dissipation and 
enstrophy induced by such vortices are equal in the limit of infinite Reynolds number. 



Recent direct numerical simulations (DNS) [y have re- 
vealed spatial structures with regions of intense dissipa- 
tion density e = li'{dui/dxj + duj/dxi)"^ and enstrophy 



density ft 



where i/ is kinematic viscosity, a; = V x u 



is vorticity, and u is the fluid velocity. The spatial distri- 
butions of enstrophy density and dissipation density are 
qualitatively different. Intense vorticity has a tube- like 
(or filament-like) structure, while intense dissipation typ- 
ically surrounds the vortex tubes and forms double-peak 
structures centered on the tubes. This difference of struc- 
ture implies differences in intermittency between enstro- 
phy and dissipation and has been linked to a difference 
between the empirical scaling exponents of longitudinal 
and transverse structure functions recently observed in 
experimental measurements and DNS 

Siggia 0, Kerr and Meneveau et al. have 

noted that enstrophy is more intermittent than dissipa- 
tion. More recently, Chen et al. |^ used DNS data to 
calculate empirical scalings of enstrophy and dissipation. 
They found that the exponents for locally-averaged en- 
strophy were significantly smaller than those for locally- 
averaged dissipation and related the differences to an 
elaboration of Kolmogorov's refined similarity hypoth- 
esis. 

Traditional experimental measurements typically in- 
voke Taylor's hypothesis and the calculation of inertial- 
range scaling exponents can be contaminated by large- 
scale shear effects. DNS using the largest available com- 
puters can reach only moderate Reynolds, at which the 
inertial range is narrow. These difficulties accentuate the 
need for theoretical understanding of intermittency phe- 
nomena and their role in scaling. 

Distributions of idealized vortices have been used by a 
number of authors to model turbulence statistics. The 
hope here is that the vortex distribution can character- 
ize both the essential physics of the small-scale structures 
and the observed vorticity statistics, thereby leading to 
an acceptable description of overall turbulence proper- 
ties. Idealized vortices that have been studied include the 
Hill vortex , simple vortex filaments |Q , the Burg- 
ers vortex ]14| , and spiral vortices [Q. A distribution 



of modified Burgers vortices will be used in the present 
paper to model scaling exponents over finite ranges. Be- 
fore presenting detailed results for the Burgers vortex, we 
shall discuss the underlying question of whether cylindri- 
cal vortex structures can support different asymptotic 
scaling exponents for enstrophy and dissipation. 

Moments of locally averaged enstropy and dissipation 
densities in a velocity field with isotropic overall statistics 
may be defined by 



(1) 



where ( ) denotes ensemble average over the isotropic 
statistics and [ ]i denotes space average over a region of 
characteristic linear dimension £. Joint powerlaw scaling 
with enstrophy exponents Cn ^^nd dissipation exponents 
^„ exists if there is a range of £ in which 



r!„(^)K(i/^)^", e„(^) (X (L/£)^ 



(2) 



Here L is a macroscale that marks the bottom of the 
scaling range. The space averaging [ ]i smears spots of 
intense excitation. Therefore C„ and ^„ are expected to 
be positive. Let £d mark the top of the scaling range, 
which we assume is also the beginning of the dissipation 
range. 

Suppose that < ^„ < Cn and that the length of the 
scaling range in decades becomes infinite in the limit 
Re — > OO where Re is a turbulence Reynolds number. 
Then f2n(^d)/en(^ci) —> oo as Re — > oo. A corollary 
is that the ratio of single-point averages becomes infi- 
nite in the limit: (fi") / (e") oo. We assume here 
that (e") jtnijld) does not become infinitely larger than 
(fi") /r2„(£d)- Now we can ask what kind of vortex struc- 
ture, if any, can support this behavior. If (il") / (e") goes 
infinite, where the averages are over the entire field, then 
there must be at least one vortex structure in the field 
for which the ratio of averages over that single structure 
goes infinite. 

A cylindrical vortex is characterized by the azimuthal 
velocity vsir), where r is distance from the vortex axis. 
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The enstrophy density and dissipation density (normal- 
ized by v) are 

The enstrophy and dissipation per unit length of vortex 
are 27r f2(r)r dr and 27r e(r)r dr, respectively. The 
moments 



17„ = 27r/ [n{r)]"rdr, e„ = 27r 
Jo 



i{r)]"rdr (4) 



describe the distribution of enstrophy and dissipation 
densities in the single vortex structure. By @, fii = ei 
if veiO) = vei^) = 

Consider first the Rankin vortex of radius rp for which 
ve oc r {r < tq) and ve oc 1/r (r > tq). Here the vor- 
ticity is confined to a rigidly rotating core and all the 
dissipation lies outside the core. The ratios i?„ — iln/e„ 
calculated from (||) are i?„ = 2n — 1. 

For the Burgers vortex of radius rj,, 



Voir) 



r 1 



exp(. 



2tt 



(5) 



where F is the total circulation. By (||), ri(r) oc 
exp(— 2r^/r^), and the vorticity and dissipation now 
overlap. Eq. (|) yields R2 ~ 10.65, R3 « 104.07, 
R4 « 1040.02. To a rough approximation, i?„ ^ 10""^. 

The cylindrical vorticity distribution that maximizes 
i?2 can be found by solving the associated variational 
problem, with u)(r)r dr and Oi held constant. The 
result confirms what can be guessed by inspection of (||): 
The maximizing distribution is the limit rg — > 00 of 

ve{r) (X r [r < rg), VQ{r) oc r" {tq < r < ri), 



ve{r) cx 1/r (r > ri). 



(6) 



with a = 1/2. In this case, i?„ = (n-l)9" for n > 1. The 
essential facts here are first, if ve oc r", then f7(r)/e(r) 
grows as a ^ 1; second, a = 1/2 is the largest a for 
which J^^ [il,{r)]^r dr diverges as ro — > 0, thereby making 
the relative contribution to 62 from the vg oc 1/r region 
negligible in the limit. The maximum R„ for n > 2 also 
are finite; they are maximized at a values that increase 
with n. 

It follows that differing asymptotic scaling exponents 
for enstrophy and dissipation cannot be supported by 
cylindrical vortices. The vortex profile that maximizes 
i?2 is one in which w(r) is highly diffuse and the domi- 
nant contributions to [ri(r)]^ and [e(r)]'^ have the same r 
dependence. This refutes the intuition that the maximiz- 
ing distribution is one in which a compact vortex core is 
surrounded by diffuse dissipation. 



Similar arguments establish that differing asymptotic 
scaling exponents for longitudinal and transverse struc- 
ture functions |^ cannot be supported by cylindrical 
vortices. 

If the cylindrical vortex is replaced by a plane vortex 
layer, it is trivial that all i?„ = 1. This suggests that the 
cylindrical vortex may be the form that maximizes the 
Rn- In any event, we conjecture that the i?„ are finite 
whatever the shape of the finite volume that is filled with 
the intense vorticity. 

The Navier-Stokes (NS) equation so far has not been 
invoked. Highly diffuse cylindrical vortices are not ex- 
pected to survive under the NS equation; the value 
i?2 = 10.65 for the Burgers vortex plausibly is closer to 
the maximum R2 attainable under NS. 

In the remainder of this paper, a distribution of mod- 
ified Burgers vortices will be constructed in order to 
demonstrate that different finite-range scaling exponents 
for enstrophy and dissipation can be realized despite the 
fact that the infinite-range exponents must be equal. 
Data from experiments and DNS are invoked in the con- 
struction, but the model is essentially artificial and we 
do not offer it as a realistic representation of the vortex 
structure of actual turbulence. 

Exponential-like distributions of vortex core radii have 
been observed in DNS by Jimenez et al. [ p^ and in an 
experiment by Tabeling et al. |^l[]. We introduce a dis- 
tribution of core sizes into our model by fitting the high- 
Reynolds number data in pl| to a stretched exponential, 
thereby obtaining the radius pdf 



P{rb) oc r] ^rl exp[ 



0.71 



(7) 



where 77 is an assumed Kolmogorov dissipation length 
p2| . Of course the experimentally observed vortex cores 
are not Burgers vortices. 

The Burgers vortex (|^) is in equilibrium under the NS 
equation if it lies within a uniform background strain field 
with radial and axial strain components —a and a, such 
that rt = (2zy/a)i/2_ jjere v is kinematic viscosity. This 
background field has no vorticity but it does contribute 
to the total dissipation density, which becomes 



dvg 
dr 



r 



(8) 



We consider a distribution of isolated Burgers vortices, 
each in its own background field. The ensemble average 
( ) is over the radius pdf (^. The local averaging [ ]i 
invoked in is realized as the average over a cylinder 
of radius £ and length 2£, centered on on a point x in 
the space surrounding a vortex and aligned with the axis 
of the vortex, followed by an average over x. In order to 
have finite total dissipation per unit core length according 
to (^) , the integration over x is extended only to a finite 
distance rmax from the vortex axis. This truncation is in 
lieu of taking account of vortex interactions and making 
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realistic changes in the straining field far from the vortex 
core. 
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FIG. 1. Parametric plot of en{£) as a function of e2{£) for 
n = 3,4,5 and 6 (from top to bottom). Finite-range ex- 
ponents used in Figs. 2 and 5 are best fits over the region 
between the two dashed lines. The inset shows e2{£) versus £. 



In the numerical study reported here, the space inte- 
grations needed to evaluate the n„ (l) and e„ {£) were per- 
formed by a second-order trapezoidal scheme whose ac- 
curacy was verified. We have taken rmax = lOry and have 
fixed the vortex Reynolds number at Rr = T /v = 1293. 
These values, together with (non- uniquely) make the 
inhomogeneous model satisfy the homogeneity relations 
= ei(^). By (||), a change in the value of i?r does 
not affect the enstrophy scaling, but by (||) it does affect 
the relative scalings of enstrophy and dissipation. 
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FIG. 2. The scaling exponents ^„ of the locally averaged dis- 
sipation as a function of n: The circles represent our Burgers 
vortex model, the diamonds are data from an experiment by 
Sreenivasan et al. and the crosses are DNS data 

Fig. 1 shows the e„(£) plotted against e2(^) for n = 
3, 4, 5. A narrow powerlaw range can be identified. In 



the inset, we plot e2(£) against £. Local scaling ex- 
ponents for dissipation and enstrophy are calculated as 
in = din e„(^)/dln £ and Cn = din VLn{l) / d\n i. The 
data used for Fig. 2 yield 6 = 0.195 ± 0.0747, which is 
not far from the value 2/9 given by some simple scaling 
models or the experimental value 0.25 ± 0.05 p3[ |. 

Fig. 2 shows in versus n. For n < 2, the present model 
agrees well with the DNS and experimental data. For 
n > 2, the exponents from the present model decrease 
too slowly, indicating an unrealistic rate of increase of 
intermittency with decrease of scale. 
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3. Parametric plot of Q,2{£) as a function of ^2{t)- 



The circles denote the part of the plot for which I lies be- 
tween the dashed lines of Fig. 1. The dotted line corresponds 
to ^2{(.) ~ ^2((-)- The inset is a similar plot of 0,4,(1) and 




FIG. 4. Ratio of enstrophy and dissipation scaling exponents 
Cn/£,n as a function of n. The circles represent the present 
Burgers' vortex model, the diamonds are data from an exper- 
iment by Antonia et al and the crosses are DNS data Q. 



In Fig. 3, fl2{£) and il4{£) are plotted against e2{£) 
and e4{£), respectively. Fig. 4 shows the ratio Cn/£,n as a 
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function of n for 4/3 < n < 4. Data are presented for ex- 
periment 1^, DNS §], and the present model. The DNS 
was computed at mesh resolution 512'^. The DNS ratio 
values in Fig. 4 decrease very slowly with increase of n. 
For the experimental data, the ratio values were calcu- 
lated from measurements of longitudinal and transverse 
velocity structure functions by invoking refined similarity 
hypotheses for longitudinal velocity increments [ p4[ and 
transverse velocity increments Experiment, DNS, 
and model all yield Cn/6i > 1- 

We reach two principal conclusions. First, if the 
support of intense vorticity lies in cylindrical vortices, 
then differing asymptotic scaling exponents for locally- 
averaged enstrophy and dissipation are impossible. We 
do not see how this conclusion can be changed if the sup- 
port lies in vortex structures of non-cylindrical shape. 
Second, models built on cylindrical vortices can yield dif- 
fering finite-range scaling exponents |2^. It is plausible 
that compact cylindrical vortices do mediate the entro- 
phy intermittency measured at small scales. Moffatt et 
al found, from an analysis of cylindrical vortices at 
a uniformly distributed angle to a constant strain field, 
that the vortex cores contained 63.3% of the total enstro- 
phy and only 1.3% of the total enstrophy. 

The Burgers vortex model presented here is arbitrary 
in a number of respects. Nevertheless, it is plausible that 
the filamentary vortex cores of actual turbulence are asso- 
ciated with moment ratios i?„ that have orders of magni- 
tude similar to those of the Burgers vortex. In this event, 
different best-fit scaling exponents for enstrophy and dis- 
sipation can be induced over substantial finite ranges of 
scales. 

We thank H. Chen, M. Nelkin, K. R. Sreenivasan, Ste- 
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